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Abstract. 

By definition, brown dwarfs never reach the main-sequence, cooling and dimming over their 
entire lifetime, thus making substellar models challenging to test because of the strong depen- 
dence on age. Currently, most brown dwarfs with independently determined ages are companions 
to nearby stars, so stellar ages are at the heart of the effort to test substellar models. However, 
these models are only fully constrained if both the mass and age are known. We have used 
the Keck adaptive optics system to monitor the orbit of HD 130948BC, a brown dwarf bi- 
nary that is a companion to the young solar analog HD 130948A. The total dynamical mass of 
0.109±0.003 Mq is the most precise mass measurement (3%) for any brown dwarf binary to date 
and shows that both components are substellar for any plausible mass ratio. The ensemble of 
available age indicators from the primary star suggests an age comparable to the Hyades, with 
the most precise age being 0.79^15 Gyr based on gyrochronology. Therefore, HD 130948BC 
is unique among field L and T dwarfs as it possesses a well-determined mass, luminosity, and 
age. Our results indicate that substellar evolutionary models may underpredict the luminosity 
of brown dwarfs by as much as a factor of «2-3x . The implications of such a systematic error in 
evolutionary models would be far-reaching, for example, affecting determinations of the initial 
mass function and predictions of the radii of extrasolar gas-giant planets. This result is largely 
based on the reliability of stellar age estimates, and the case study of HD 130948A highlights 
the difficulties in determining the age of an arbitrary field star, even with the most up-to-date 
chromospheric activity and gyrochronology relations. In order to better assess the potential 
systematic errors present in substellar models, more refined age estimates for HD 130948A and 
other stars with binary brown dwarf companions (e.g., e Ind Bab) are critically needed. 
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1. Introduction 

Theoretical models of objects below the substellar limit (M < 0.075 M Q ) are essen- 
tial for characterizing the several hundred brown dwarfs and extrasolar gas-giant planets 
discovered to date. Thus, these models have become ubiquitous in the literature, even 
though empirical tests of their ability to accurately predict the properties of brown dwarfs 
has been limited to only a handful of relatively warm objects. To test substellar evolu- 
tionary models, the input parameters of mass and age must be determined. For young 
brown dwarfs, the M6.5 eclipsin g binary 2MASS J05 352184-0546085 in the Orion Nebula 
provides a unique benchmark ( Stassun et al. 200fiJ). Prior to this year, only three bina- 
ries provided dynamical mass measuremen ts for field objects at or below the sub stellar 
limit: the M8.5+M9 binary L HS 1070BC dLeinert et aljfeoOlt ISeifahrt et all [20081) : the 
M8.5+M9 binary Gl 569Bab (jZapatero Osorio et all 12004 ISimon et all 120061) : and the 



2 Dupuy et al. 

L0.5+L1 binary 2MASS J0746+2000AB (jBouv et alj|2004) . Recent work has contributed 
several more dynamical masses for objects lo wer in both temper ature and mass than 
previously studied: the mid-L dwarf GJ 802B (llreland et al.|[2008h : the T5+T5.5 dwarf 
binary 2MASS J15 34-2952AB (|Liu et al.ll2008h ; and the L4+L4 binary HD 130948BC 



( Dupuv et al.l 120081 ). While mass measurements alone can provide very stringent tests 
of theoretical models (e.g., see Liu et al. 2008), substellar evolutionary models are only 
fully constrained when both the mass and age can be determined. In fact, precise ages are 
critical for such tests because brown dwarfs - unlike stars - never reach a main-sequence, 
so their properties depend very sensitively on their age. 

Of the substellar field dwarfs with measured masses, only HD 130948B C has a pre- 
cisely determined age. These nearly-identical L dwarfs were discovered by IPotter et al.1 
(|2002l ) as companions to the young solar analog HP 130948A (G2V [Fe/H] = 0.05). 
Hipparcos measured a distance of 18.17±0.11 pc ( van Leeuwen 2007f ) for the primary 
star, which enables a very precise dynamical mass measurement when paired with our 
wcll-dctcrmincd orbital solution. 



2. The Mass of HD 130948BC 

We have used Keck adaptive optics (AO) imaging to monitor the relative orbit of 
the two components of HD 130948BC (Figure [JJ). Combined with archival Hubble Space 
Telescope (HST) imaging and a re-analysis of the Gemini discovery data, our data span 
«7 years («70% of the orbital period). We fit a simple analytic PSF model to derive 
astrometry from the Keck and Gemini images, while TinyTim model PSFs were fit to the 
HST images. An individually tailored Monte Carlo simulation was used to determine the 
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Figure 1. Keck (circles), HST (squares), and Gemini (triangle) relative astrometry for 
HD 130948BC along with the best-fit orbit. Error bars are smaller than the plotting symbols. 
The empty circles are the predicted positions in 2009 and 2010. 
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astromctric uncertainty for each observation epoch. The resulting astrometry is extremely 
precise with typical Keck errors of 300 fias, corresponding to ml R Q at the distance of 
this system, while the orbit is roughly the size of the asteroid belt. We determined the 
binary's orbital parameters and their confidence limits using a Markov Chain Monte 
Carlo (MCMC) technique. The best-fit orbit has a reduced x 2 °f 1-06 (9 degrees-of- 
freedom), thus validating our astrometric error estimates. Applying Kepler's Law to the 
MCMC-derived orbital period (P = 9.9±° r 6 yr) and semimajor axis (a = 121 ± 6 mas) 
yields a dynamical mass of 0.10891q'qqi7 Mq. Accounting for the additional uncertainty 
in the Hipparcos distance results in a dynamical mass of 0.109±0.003 M Q (114±3 Mj up ). 

In the following analysis, we apportion the total mass between the two components by 
converting the measured luminosity ratio into a mass ratio using evolutionary models. 
The resulting individual masses are very insensitive to the models used because the flux 
ratio is so close to unity (the steepness of the mass-luminosity relation means that even 
small differences in mass result in large differences in luminosity). Regardless, we are 
careful to conduct our analysis in a self-consistent manner free of circular logic. 



3. The Age of HD 130948A 

As a young solar analog, multiple indicators are available to assess the age of HD 130948A: 



• Rotation/ Gyrochronology — Gaidos et al. ( 2000h measured two rotation periods of 
7.69 and 7.99 days for HD 130948A. Thus, we adopt a rotation period of 7.84±0.21 days 
and a B — V color of . 576±0 .016 mag from the Hipparcos catalog. We employ the 
Mamaiek fc Hillenbrand J2 008) calibration of the "gyrochronology" relation originally 



introduced by iBarnesI (|2007l ). The age we derive is 0.79_q 15 Gyr, where the confidence 
limits are determined through a Monte Carlo approach in which the period, color, and 
empirical coefficients are drawn from normal distributions corresponding to their uncer- 
tainties. 



Henry et all (|l996h and lWright et all (|2004l ) measure log(i?; 



• Chromo spheric Activity 
v alues of —4.45 and —4.50 f or HD 130948A, respectively. Using the activity-age relation 
of lMamaiek fc Hillenbrand! (|2008h . we derive ages of 0.4 and 0.6 Gyr from these \og(R' HK ) 
values. The empirical relation is expected to gives ages with an uncertainty of «0.25 dex, 
so we adopt a mean age of 0.5±0.3 Gyr from this method. 



hk ) 



• X-ray Activity 
measure \og(Lx) 



(1999) 



- HD 130948 A was detected by ROSAT, and iHiinsch et al. 
29.0 1 dex (cgs), which gives log(i?x) = —4.70. Using the empirical 



relation of lGaid os (1998), this corresponds to an age of 0.1-0.3 Gyr, depending on w hether 
we adopt a of 0.5 or 1/ exp. The X-ray relation of lMamaiek fc Hillenbrand ( 2008 ). de- 
rived by combining their log(Rx)-^og(R' H K ) and \og(R' HK )-age relations, gives an age 
of 0.5 Gyr. The X-ray luminosity of HP 130948A is in agreement with single G star s in 
the Pleiades and Hyades (28.9-29.0; fstern et all (|l995l ); IStelzer fc Neuhauserl <|200lh V 



• Isochrones — Using high resol ution spectroscopic data combined with a bolometric 
luminosity and model isochrones, IValenti fc Fischer! (|2005h derived an age estimate of 
1.8 Gyr, with a p ossible age range of 0.4-3.2 Gyr. From the same data and with more 



detailed analysis, Takeda et al.1 ( 2007 ) found a median age of 0.72 Gyr, with a 95% 



confidence range of 0.32-2.48 Gyr. 

• Lithium — Measurements bv lDuncanl ()l98lh . lHobbj (|l985h . and lChen et all (|200lh give 



lithium equivalent widths of 95±14, 96±3, and 103±3 mA, respectively, for HD 130948A. 
Compared to stars of similar color, these values are slightly lower than the mean for the 
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Figure 2. The filled circles mark the measured luminosities of HD 130948B and C at the age 
we derive for HD 130948A. The thick shaded lines are isomass lines from evolutionary models, 
where the line thicknesses encompass the la errors in the individual masses of HD 130948BC. 
Although the two independent sets of models agree very well with one another, they underpredict 
the luminosities of HD 130948BC by a factor of «2-3x. 



Pleiades and slightly higher than for UMa and the Hyades, tho ugh consistent with the 
scatter in each cluster's measurements ISoderblom et al. (|l993al fbllcl) 



In summary, the most precise age estimate available for HD 130948A comes from 
gyrochronology, which gives an age of 0.79^Q^g Gyr. All other age indicators agree with 
this estimate, though this is due to their large uncertainties rather than a true consensus. 



4. Substellar Evolutionary Models Fully Constrained 

With a measured mass, luminosity, and age, HD 130948BC provides the first direct 
test of the luminosity evo lution predicted by t heoretical models for substel lar field dwarfs. 
Both the Tucson models (I Burrows et al 1ll997t) and Lyon models (DUSTY; [Chabr ier et al 



2000) underpredict the luminosities of HD 130948B and C given their masses and age. 
The discrepancy is quite large, about a factor of 2 for the Lyon models and a factor of 3 for 
the Tucson models (Figure [2]). If the age and luminosities of HD 130948B and C had been 
used to infer their masses, the resulting estimates would have been too large by 20-30%. 
In order to explain this discrepancy entirely, model radii would have to be underpredicted 
by 30-40%. Alternatively, the age of HD 130948A would need to be «0.4 Gyr in order 
to resolve this discrepancy. Although such a young age is marginally consistent with 
the various age indicators; it is on the extreme young end of two independent, well- 
calibrated age estimates (gyrochronology and stellar isochrones). In order to better assess 
this discrepancy between models and data, a more refined age estimate for HD 130948 A 
(e.g., from asteroseismology) is critically needed. 



5. Lithium Depletion in HD 130948BC 

Since brown dwarfs are fully convective objects, they can rapidly deplete their initial 
lithium if their core temperature is ever high enough to do so. This threshhold is reached 
around 0.065 M Q , and since this is below the hydrogen-burning mass- limit, this fact has 
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Figure 3. Lithium depletion as a function of age as predicted by evolutionary models. The solid 
lines correspond to the individual masses of HD 130948B and C. These lines are bracketed by 
dotted lines that correspond to the ler uncertainties in the individual masses. The ordinate is 
the fraction of initial lithium remaining. The hatched black box indicates the constraint from 
the age of HD 130948A estimated from gyro chronology. 

been exploited to identify sufficiently old objects bearing lithium as substellar. In fact, 
the exact mass-limit for lithium burning is slightly different depending on which sets of 
theoretical models are used, and the masses of HD 130948B and C happen to be very close 
to these theoretically predicted mass-limits (Figure [3]). According to the Tucson models, 
neither component is massive enough to have ever depleted a significant amount of its 
initial lithium. The Lyon models, on the other hand, predict that HD 130948B is massive 
enough to have depleted most of its lithium, while HD 1 30948 C is not. Thus, resolved 
optical spectroscopy designed to detect the lithium doublet at 6708 A would provide a 
very discriminating test of substellar evolutionary models, which are otherwise nearly 
indistinguishable (e.g., see Figure El). This experiment can currently only be conducted 
with HST/STIS given the very small binary separation (< 130 mas). 

6. Future Prospects 

Brown dwarfs hold the potential to address many astrophysical problems. For example, 
they are excellent laboratories in which to study ultracool atmospheres under a variety 
of conditions, and they may eventually be useful as Galactic chronometers given how 
sensitively their properties depend on their age (see contribution by A. Burgasser). How- 
ever, the theoretical models we rely upon to characterize brown dwarfs have only begun 
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to be rigorously tested by benchmark systems such as HD 130948BC. More results are 
expected to be forthcoming over the next several years for oth er brown d warf binaries 
with stellar companion s: e Ind Bab ( McCaughrean et al. 2004 ); Gl 417BC ( Bouv et al 



2003j) ; and GJ 1001BC (jGolimowski et al.ll2004) . However, the utility of these systems as 



benchmarks critically depends on the confidence in the age estimates for their primary 
stars. Therefore, these stars deserve special attention so that state-of-the-art age-dating 
techniques (e.g., asteroseismology and gyrochronology) may be applied to them. Also, 
extending the empirical relations between age, stellar rotation, and chromospheric ac- 
tivity to include objects with as late a spectral type as possible will enable many more 
systems to be used as benchmarks for testing models. These relations are currently only 
calibrated for stars as late as early-K, but about half of the stars with brown dwarf 
companions have spectral types between early-K and early-M. 
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Discussion 

F. Walter: It is always risky to attempt to pin down the age of a field star, even using 
multiple techniques that may agree. How much would the age have to be changed to 
place the L dwarfs on the proper evolutionary tracks? 

T. Dupuy: I agree and would really like to see another independent measurement of the 
age, for example, from asteroseismology. The age of the system would have to be about 
0.4 Gyr to bring the models into agreement with the data. 

E. Jensen: Is there a measured metallicity for HD 130948A? The metallicity will affect 
the evolutionary models, both in the H-R diagram and the predicted Li depletion. 

T. Dupuy: That's exactly right; a detail I didn't go into. The metallicity of HD 130948A 
is basically solar, which means we can use the standard models. This is another reason 
why having brown dwarfs with stellar companions is great - because you can make sure 
you're not being confused by metallicity effects like those Adam talked about. 

A. West: Does the fact that this system is a close binary affect the measured luminosity 
(because the radii are affected)? 

T. Dupuy: The binary separation is about 2.2 AU, so it's unlikely that tidal effects 
are at work in this system. Also, it turns out that the two components receive about as 
much flux from each other as they do from the primary star, so irradiation shouldn't be 
affecting them much. 

J. Fernandez: The next main source of benchmarks for brown dwarfs will be Kepler 
and Corot. The precise determination of ages for these primary stars will be crucial. 



